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Abstract: We introduce a quantum stochastic dynamics for heat conduction. A multi- 
level subsystem is coupled to reservoirs at different temperatures. Energy quanta are 
detected in the reservoirs allowing the study of steady state fluctuations of the entropy 
^ . dissipation. Our main result states a symmetry in its large deviation rate function. 

KEY WORDS: entropy production, fluctuation theorem, quantum stochastic 
calculus 

C/3 



1 Introduction 



'^ ' Steady state statistical mechanics wants to construct and to characterize the 

Q ■ stationary distribution of a subsystem in contact with several reservoirs. By 

, ^, . nature the required scenario is an idealization as some essential specifications 

of the reservoirs must be kept constant. For example, intensive quantities 
such as temperature or (electro-) chemical potential of the different reservoirs 



"Tj- ' are defined and unchanged for an extensive amount of time, ideally ad in- 



finitum. Reservoirs do not interact directly with each other but only via the 



v^ . subsystem; they remain at their same spatial location and can be identified 

O i at all times. That does not mean that nothing happens to the reservoirs; 

Q I fiows of energy or matter reach them and they are like sinks and sources of 

currents that fiow through the subsystem. Concrete realizations and models 

of steady states vary widely depending on the type of substances and on the 

nature of the driving mechanism. 
'^ ! An old and standard problem takes the subsystem as a solid in contact at 

Q I its ends with two heat reservoirs and wants to investigate properties of the 

J-^ ' energy fiow. Beloved by many is a classical model consisting of a chain or an 

array of coupled anharmonic oscillators connected to thermal noises at the 
^ . boundaries. The reservoirs are there effectively modeled by Langevin forces 

while the bulk of the subsystem undergoes a Hamiltonian dynamics, see e.g. 

|H,Bl'n21 IEPRB99llMNVn3j . Our model to be specified below is a quantum 
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analogue of that scenario in the sense that we also consider a combination of 
Hamiltonian dynamics and Markovian thermal noises. 

We imagine a chain of coupled two (or multi)-level systems. The dynamics 
of the isolated subsystem is unitary with Hamiltonian Hg. Quanta of energy 
uj are associated to the elementary transitions between energy levels. Two 
physical reservoirs at inverse temperatures (3k, k = 1,2, are now attached 
to the subsystem. The total dynamics is described by a quantum stochastic 
differential equation through which we can observe the number N^^^ of quanta 
with energy u that are piled up in the k—th reservoir. The total energy 
N := Hs + "^uj k^^'^,k is conserved under the dynamics (Proposition 12. IHJ) . 
The change in the second term corresponds to the flow of energy quanta in 
and out of the reservoirs and specifies the dissipated heat. Our main result 
consists in obtaining a symmetry in the fluctuations of that dissipated heat 
that extends the so called steady state fluctuation theorem for the entropy 
production to a quantum regime (Proposition 12. l4jl . 

The quantum stochastic evolution that defines the model is a particular 
dilation of a semigroup dynamics that describes the weak coupling regime 
of our subsystem coupled to quasi-free boson fields. The dilation, a sort of 
quantum Langevin equation, is much richer and enables the introduction of a 
natural path space measure. One should remember here that a major concep- 
tual difficulty in coming to terms with the notion of a variable entropy pro- 
duction for quantum steady states is to understand its path-dependence. One 
option is to interrupt the unitary dynamics with collapses, see e.g. |CRJ"'"04J . 
Others have proposed an entropy production operator, avoiding the problem 
of path-dependence. Our set-up follows a procedure that is well-known in 
quantum optics with thermal noises formally replacing photon detectors, see 
|BGM04l [BMK03J . In the resulting picture we record each energy quantum 
that is transferred between subsystem and reservoirs. It induces a stochastic 
process on quanta transferrals and there remains no problem to interpret 
the fluctuations of the entropy production. From the mathematical point of 
view, the model can be analyzed via standard probabilistic techniques. 

1.1 Related results 

In the past decade, a lot of interest has been going to the Gallavotti-Cohen 
fluctuation relation |(;(^95bl IHTW^^ IHX^M93] . see |Ma,en3J for more recent 
references. In its simplest form that relation states that the steady state 
probabihty (Prob) of observing a total entropy decrease wt = —wT in a 
time T, is exponentially damped with respect to the probability of observing 



an increase of wt as 



Prob{wT = wT) _^+^^ 



Prob{wT = —wT) 

at least for very large time spans T. The relation (jl.lll is known as the 
steady state fluctuation theorem (SSFT) and states a symmetry in the fluc- 
tuations of the entropy dissipation in a stationary nonequilibrium state. The 
symmetry was first discovered in the context of dynamical systems and was 
applied to the phase space contraction rate in strongly chaotic dynamical 
systems, see |R,ue99| IE(yM93[ 'GC95b' 'GC95aj . It was first further developed 
for stochastic dynamics in [Kur98. .1899. ■Mae99j . In the present paper, we 
deal with the SSFT for a quantum system. The steady state condition must 
however be understood in a physical sense; it is about heat conduction for 
fixed reservoirs in the long time limit. The small system is treated in the 
steady state of an approximate dynamics (the weak coupling limit) while the 
reservoirs are kept at a fixed termperature. Yet, mathematically, we are not 
quite deahng with a steady or stationary state. The true total dynamics for 
system plus reservoirs is in fact much more complicated. We still speak about 
the steady state (and the SSFT) also to contrast it with transient versions 
of the fluctuation symmetry (jl.lj) . see also |(T(y99j . Transient fluctuation 
theorems (TFT) start typically from a change of variables at a finite time t, 
reversing so to say the evolution, see JSEC98J . and can be obtained equally 
well for classical as for quantum systems. That is not at all what we are 
doing here. We truly concentrate on the stationary heat dissipation in the 
reservoirs, but from a technical point of view, one could argue that our set-up 
is actually a "transient model of a steady state" . 

The basic underlying mechanism and general unifying principles connect- 
ing SSFT and TFT with statistical mechanical entropy have been explained 
in pNfM lM^^. 

Monnai and Tasaki |MTn8bj have investigated an exactly solvable har- 
monic system and found quantum corrections to both SSFT and TFT. Matsui 
and Tasaki |MT03aj prove a quantum TFT in a general C*-algebraic setting. 
It is however unclear what is the meaning of their entropy production oper- 
ator. 

A related quantum Jarzynski relation was studied in |RMn4j . 

Besides the fluctuation theorem, we also describe a new approach to the 
study of heat conduction in the quantum weak coupling limit. In |LS78j . 
Lebowitz and Spohn studied the thermodynamics of the weak-coupling gen- 
erator. They identified the mean currents, and they proved a Green-Kubo 
relation. At that time it was however not yet possible to conclude that these 
expressions are the first non-zero contributions to their counterparts at finite 



coupling A. That has recently been shown in a series of papers by Jaksic and 
Pillet [JP96a . IJPQHH 1,1 Pn2 bl lJPn 2aj . who used spectral techniques to study 
the system at finite coupling A. It was also shown that the stationary state 
of the weak coupling generator is the zeroth order contribution to the system 
part of the so called NESS, the nonequilibrium steady state. The current 
fiuctuations we define in our model, agree with the expressions of |LS78j as 



far as the mean currents and the Green-Kubo formula is concerned. Our 
entropy production operator is however new; it differs for example from the 
proposal of |MTn8aj . The approach taken here also differs from the more 
standard route that has been followed and that was outlined by Ruelle in 
jRueOOj . Recently and within that approach and context of heat conduction, 
new results have been obtained in |JP02b| IJP02al IASF05J . To us it remains 
however very much unclear how to define and study in that scenario a fiuc- 
tuating entropy; in contrast, that is exactly one of the things we can easily 
achieve via our approach but we remain in the weak coupling limit. 

1.2 Basic strategies 

1.2.1 Microscopic approach 

In general one would like to start from a microscopic quantum dynamics. 
The system is then represented by a finite-dimensional Hilbert space Ti and 
system Hamiltonian Ho,. The environment is made from thermal reservoirs, 
indexed by /c G -ft', infinitely extended quantum mechanical systems, with 
formal Hamiltonian, 

keK 

The coupling between system and reservoirs is local and via some bounded 
interaction term Ai^s-R so that the total Hamiltonian takes the form 

keK 

where we have already inserted a specific form for the coupling -f^s-R using 
self-adjoint reservoir operators Rk and Vk acting on respectively TYr^. and H. 
On the same formal level, which can however easily be made precise, the 
total quantum dynamics is then just 



Ut 



A .^ ^~iH),t 



We will not follow the beautiful spectral or scattering approach that has 
recently been exploited for that nonequilibrium problem. We refer the reader 



to the specialized references such as |JP02bl IRueOOj and we only outline the 
main steps, totally ignoring essential assumptions and technicalities: 
One starts the dynamics from an initial state 



At :=Ps ®PRi ® • •• ®PRj 



A' I 



where ps stands for an initial state in the system and the pr^ are equilibrium 
KMS states at inverse temperature j3k for the fc-th reservoir. The quantum 
dynamics takes that initial state to the new (now coupled) state pj at time 
t > 0. The NESS is obtained via an ergodic average 



1 '■^ 



Pness := lim — / dt ^t (1-2) 

TT+oo T Jo 

One of the first questions (and partially solved elsewhere, see e.g. |,TPn2b[ 
IRueOOtFASOGj ) is then to derive the natural conditions under which the mean 
entropy production rate 

m 

S := iy ^(3k Pness([-^a, -^rJ) 
fc=i 

is strictly positive. While that mean entropy production certainly coincides 
with conventional wisdom, we do not however believe that the operator 

or equivalent expressions, is the physically correct candidate for the study 
of current fluctuations which would obey the SSFT. That is not even the 
case for the simplest (classical) stochastic dynamics; one needs to go to path 
space and study current fluctuations in terms of (fluctuating) trajectories. 

1.2.2 Weak coupling approach 

Starting from the microscopic dynamics above, we can of course always look 
at the reduced dynamics A^ on the system 

A^PS = Ttr [f/,\ps ® PR, ® . . . ® PR|K|)t^^i 

for a density matrix ps on the system. Obviously, the microscopic evolution 
couples the system with the environment and the product form of the state 
will in general not be preserved. One can however attempt a Boltzmann-type 
Ansatz or projection technique to enforce a repeated randomization. That 
can be made rigorous in the so called weak coupling limit. For that, one 



needs the interaction picture and one keeps \^t = r fixed. That is the Van 
Hove-Davies-hmit |Hov55| IDav74j 

limA°_,A,Vs:=e^^*ps 

A — *U 

where C* is a hnear operator acting on density matrices for the system. The 
generator will be written out more explicitly in Section 12.11 but its dual C 
acting on B(7i) is of the form, see ()2.1()|1 . 

where the Ck can be identified with the contribution to the dissipation from 
the fc'th reservoir. Hf is an effective, renormalized Hamiltonian depending 
on details of the reservoirs and the couphng. 

From now on, we write p for the (assumed) unique invariant state (see also 
Remark EHJ: 

e^^*p = p, r > 

Again one can study here the mean entropy production, as for example done 
in |LS78j and argue that 

TApCkHs] 

represents the stationary heat flow into the fc'th reservoir, at least in the 
weak coupling regime. Nothing tells us here however about the physical fluc- 
tuations in the heat current for which higher moments should be considered. 
In fact, the reservoirs are no longer visible as the weak coupling dynamics 
is really a jump process on the energy levels of the system Hamiltonian, see 
further in Section Kill The heat flow and the energy changes in the individ- 
ual reservoirs cannot be reconstructed from the changes in the system. The 
present paper uses a new idea for the study of the fluctuations of the heat 
dissipation in a reservoir. 

1.2.3 Dilation 

While the weak coupling dynamics is very useful for problems of thermal 
relaxation (one reservoir) and for identifying the conditions of microscopic 
reversibility (detailed balance) characterizing an equilibrium dynamics, not 
sufficient information is left in the weak coupling limit to identify the variable 
heat dissipated in the various reservoirs. Heat is path-dependent and we need 
at least a notion of energy-trajectories.^ The good news is that we can obtain 

^At least, if one has a stochastic or effective description of the system dynamics, as is 
the case in the weak-coupling limit. We do not claim at all that the trajectory-picture is 
microscopically fundamental. 
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such a representation at the same time as we obtain a particular dilation of 
the weak coupling dynamics. The representation is basically achieved via an 
unraveling of the weak coupling generator L and the corresponding Dyson 
expansion of the semigroup dynamics. That will be explained in Section 1221 
There are many possible dilations of a quantum dissipation. It turns 
out that there is a dilation whose restriction to the system coincides with 
the Dyson representation in terms of energy-trajectories of the weak cou- 
pling dynamics. That dilation is well studied and goes under the name of 
a quantum stochastic dynamics. The associated quantum stochastic calcu- 
lus was invented by Hudson and Parathasaraty, |HP84j . It has been ex- 
tensively employed for the purpose of quantum counting processes, see e.g. 
|BGM04[ iBMKOSj. Various representations and simplifications have been 
added, such as in |AF87j where a (classical) Brownian motion extends the 
quantum dissipation. Unravelings of generators have been first employed in 
quantum optics in |SD81j . they are further discussed in |(yar93j . 

1.2.4 Results 

We prove a symmetry in the large deviation generating function of the dis- 
sipated heat f Proposition 12. 8p . This function is analytic and this implies 
the large deviation principle. The symmetry is recognized as the fluctuation 
theorem for the entropy production. The precise form of the fluctuation the- 
orem depends on whether the model has been derived from a reversible or an 
irreversible (e.g. because of the presence of magnetic flelds) dynamics. This 
point was clarifled in |Mae04j . By a theorem of Bryc |Bry93| , analyticity of 
the generating function implies the central limit theorem for the currents. We 
do not stress this point but it is implicitly used in deriving a Green-Kubo re- 
lation and Onsager reciprocity (Proposition 12. 1(]|) . or modiflcations of these, 
again depending on the reversibility of the original model. In all cases, the 
fluctuation symmetry helps to estabhsh strict positivity of the entropy pro- 
duction f Proposition 12. 9p . Let us stress that our main result. Proposition 
12. 8[ depends on an interpretation, as described above under Section 11.2.31 
However, the consequences of our main result. Propositions 12.91 and 12.1(11 do 
not depend on this interpretation. This will be further discussed in Section 

m 

1.2.5 Comparison with earlier results 

Technically, our fluctuation theorem is very close to the results obtained in 
|Kur98j or |LS99J. The Green-Kubo relations and Onsager reciprocity have 
been established recently in e.g. |,T()P05j for the spin-fermion model. In the 



weak-coupling limit they were discussed already in |LS78j . however there the 



authors did not distinguish between reversible and irreversible models (this is 
commented upon in Remark |2.6|) . The strict positivity in the weak-coupling 
limit was proven in |,TPn2bj ^ (for the spin-fermion model) and in |ASn6j (un- 
der general conditions). Our theorem on strict positivity is however slightly 
more general: Assuming the existence of a unique, faithful stationary state, 
we formulate a necessary and sufficient condition for strict positivity. 

1.3 Outline of the paper 

In Section 2, we introduce the quantum stochastic model and state the result. 
In Section 3 follows a discussion where the main points and novelties are 
emphasized. Proofs are postponed to Section 4. 



2 The Model 

2.1 Weak Coupling 

We briefly introduce here the weak coupling dynamics without speaking 
about its derivation, which is not relevant for the discussion here. Some 
of that was briefly addressed in Sections 11.2.1111. 2!^ and it is covered in detail 
in |l)av74j and fL^. 

Let H be a finite-dimensional Hilbert space assigned to a small subsystem, 
called system in what follows. Let H^ be a self-adjoint Hamiltonian on 7i. 
Introduce the set of Bohr frequencies 

F := {u eR\3e,e' espHs:ui = e-e'} (2.1) 

Remark that F is the set of eigenvalues of the derivation — i[iys) "]• We label 
by fc G -ft' (a finite number of) different heat reservoirs at inverse temperatures 
Pk < oo. To each reservoir k is assigned a self-adjoint operator Vk G B(H) 
and for each k E K,uj & F, we put 

V^,k= J2 UHs)Vkle'{Hs) (2.2) 

e, e' e spHs 

to = e — e' 

where le(-f^s) for e G sp(iJs) is the spectral projection on e associated to Hs- 

^Besides, from jJP02b| . it follows that the strict positivity remains true at small nonzero 
coupling, without taking the weak coupling limit. 



Fix for k E K, nonnegative functions rjk G L^ (M) and assume them to be 
Holder continuous in F C M and satisfying the condition 

r^^{x) = e~'^'=^%(-2;) > 0, x G M (2.3) 

which is related to the KMS equilibrium conditions in the reservoir k E K, 
see further under remark 12.51 Write also ioi u E F, k E K 

Skico) = hm / ^^ (2.4) 

•^iO JlR\[uj-e,uj+€] ^ — X 

which is well defined by the assumption of Holder continuity for r/^gx- From 
now on, we simply write the indices u, k for u E F,k E K. We consider the 
self-adjoint Hamiltonian 

Hf:=Y,SkHV:^kV.,k (2.5) 

satisfying by construction 

[Hf,Hs]=0 (2.6) 

We work with the following generator C on B{7i) 

£(■) = i[Hj, .] + J2vki^){v:^, ■ V.,k - l{V:,kV.,k, ■}) (2.7) 

Lij,k 

Putting T(H) C BiH) the set of all density matrices on H, i.e: 

/i e T{n) ^ Tr [/i] = 1, /i > (2.8) 

one introduces the dual generator C* on T(7Y), defined through 

Tr [AC*ij] = Tr [fiCA] , A e 0(H), /x G T{H) (2.9) 

By grouping all terms with the same k in ()2.7p . we can also write 

£(-) = -i[if/,-] + $^A(-) (2.10) 

keK 

Both C and £fcex ^ire of the Lindblad form |Lin75j and hence they generate 
completely positive semigroups e^^ and e*^''. A p E TiTi) is a stationary 
state for the semigroup e*^ iff 

£*p = or, equivalently, e^'' p = p (2-11) 



We fix an anti-unitary operator T on 7i, which has to be thought of as playing 
the role of time reversal. Let 

Hl:=THsT, V,' := TV^T (2.12) 

That defines a new model, satisfying all necessary requirements. This model 
can be thought of as the time-reversal of the original one. 

We will need the following assumptions: 
Assumption Al 

We ask triviality of the commutant 

{vl^\iu)V^,k \k e K,uje FY = CI (2.13) 

where for A C BiTi), 

Sg^'^VAg^: [A,fi] =0 (2.14) 

That ensures the existence of a unique stationary state, as stated in Remark 

Assumption A2 

We ask that the system can complete a closed cycle in which the entropy 
production is nonzero. More precisely, there are sequences a;i, . . . ,a;„ in F 
and ki, . . . ,kn in K such that 



Y^Pk^^il^O (2.15) 



i=l 



2. There is a one- dimensional projection P G B{7i) such that 

Tr [PK„,,„ . . . VL,,fe,K.„fe,P] ^ (2.16) 

Assumption A3 

This assumption expresses that our model is time-reversal invariant. It will 
be used in deriving the full fluctuation theorem, the Green-Kuho relations and 
Onsager reciprocity. 

Hi = Hs, ykeK: V^ = Vu (2.17) 
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Remark 2.1. If Assumption Al holds, then, by a theorem of Frigerio (The- 
orem 3.2 in |Fri78p and the fact that PkeK < oo, the semigroup e*^ has a 
unique stationary state p. This state is faithful, i.e. for all nonzero projections 
P^Oe Bin): 

Tr[pP]>0 (2.18) 

Assumption Al is actually a necessary condition for the existence of a unique 
stationary state. 

Remark 2.2. Assumption A2 comprises the intuitive assumption that the 
system does not break up in independent subsystems which are coupled sep- 
arately to the reservoirs. If that would be the case, then most of our results 
still hold but they become trivial. For example, the rate function e from 
Proposition 12.81 satisfies V/t G C : e(K) = 0. 

Remark 2.3. If for all k G K, Pk = P for some P, then 

pp := exp{-(3Hs)/TT [exp{-f3Hs)] (2.19) 

is a stationary state for e*^, as follows from the condition ()2.3|) and the 
explicit form ()2.7|) . 

Remark 2.4. If Al holds (assuring the uniqueness of the stationary state), 
then one easily checks 

yAeB{n),\/e^e' espHs: lim e*^(le(i/s)^le'(^s)) = (2.20) 

which is usually called "decoherence" . As a consequence of ()2.2()j) . the sta- 
tionary state p G TiTi) of e*"^ satisfies, 

Y, UHs)pUHs)=p (2.21) 

e£spHs 

Remark 2.5. If one would derive the model from a microscopic setup, then 
we can be more specific. Let Ti/j^. be the Hilbert space of the fc'th reser- 
voir and pk a thermal equilibrium state at Pk on (a subalgebra of) B(TCRf.). 
Assume the coupling is given by 

J^Vk^Rk, Rk = Rl&B{HR,) (2.22) 

k&K 

Then the functions rjk are fourier transforms of the autocorrelation function 
of Rk and the KMS conditions imply ()2.3|1 . All this is discussed at length in 
|L?99, . 

The restriction to couplings of the form ()2.22|) . where each term is self- 
adjoint by itself, is not necessary. Besides, one can also have multiple cou- 
plings per reservoir. Since this complicates our notation without introducing 
any novelty, we adhere to the simple form ()2.22j) . 

11 



Remark 2.6. If H^ is nondegenerate, one can choose T as follows: Let 
ipe, ^ ^ spi^s be a complete set of eigenvectors for Hs and put 

T{ ^ Ce^e) = 5Z ^^'^^' Cegsp/^s ^ ^^ (2-23) 

e£spHs e£spHs 

Although this does not necessarily imply Assumption A3, it does imply 

Hl = Hs, V^,,AV:^, = Vl.AiV^.y AeB{H) (2.24) 

which, as one can check from the proofs, can replace A3 for all purposes 
of this paper. Hence, a nondegenerate model is automatically time-reversal 
invariant. This explains why in |LS78j the Green-Kubo relations were derived 
for nondegenerate Hamiltonians without speaking about microscopic time- 
reversal. It also explains why time-reversal does not appear naturally in the 
framework of classical Markov jump processes. 

2.2 Unraveling the generator 

We associate to that semigroup dynamics, generated by ()2.7|) . a pathspace 
measure by a procedure which is known as "unraveling the generator" . Ba- 
sically, we will introduce \F\ x 1^1 Poissonian clocks, one for each reservoir 
and each Bohr frequency. Whenever clock (u;, k) ticks, our system will make 
a transition with Bohr frequency uj, induced by reservoir k. This will be our 
'a priori' measure da (see further). If Ho, is nondegenerate, then it is very 
easy to upgrade da to the appropriate pathspace measure: one multiplies da 
with a certain factor for each jump and with factors for the waiting times, 
obtaining something of the form 

dPp„(a) = e-(*-*")'^"+ic„. . .C2e-(*2-*^)"icie-*^''Ma (2.25) 

for some positive numbers ci, . . . , c„ and ri, . . . , r^+i and initial state po- 

When Ho, is degenerate, one has to do things more carefully, leading to 
the expression ()2.42|) in Lemma [2. 71 The technical difference between degen- 
erate and non-degenerate H^ is further discussed in Section El 



2.2.1 Preliminaries 

Put 

J7ji := {a C [0,t]||a| < oo}, fi^ := {a C M+||a| < oo} (2.26) 
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where |cr| is the cardinahty of the set cr C M. Let (fij)^^^, 1]^^ ^ stand for 
identical copies of Q], Q^ and put 

Q:= X fii „ Q, := X (fii),,, (2.27) 

O and fit are called Guichardet spaces, see |(Tui72j . An element a & fl looks 
like 

a = {ui, ki, ti; . . . ; Un, kn, tn) with < ti < t2 < • • • < ^n < +00 (2.28) 

Alternatively, and corresponding to the product in ()2.27|1 : 

(y = icroj,k)iu,k with a^^k e Vll^^, \a\ := ^ |(T^,fc| (2.29) 

LO.k 

We define integretation on ^2^ and Q, by putting for any sequence of functions 
g = {gn)nen with gn a measurable function on F^ x K^ x (R+)" for all n G N, 

/ d(T5f((T) := V] y2 / dti...dt„5f„((u;i,A;i,ti;...;u;„, A;„,t„)) 

L^i,...,w„eF" 

(2.30) 
where A" C M" is the simplex 

(ti, . . . , t„) G A^' ^ < ti < . . . < t„ < t (2.31) 

The equality ()2.30p defines the symbol "dcr" and the notion of measurable 
sets in Qt or Q (for the latter, take t = cxd in the above definitions). 
For future use, we introduce 'number functions' n^j^, defined as 

iT'lA'^) •= l^'^.fc ^ i^'^]!' '^t^.fclf^) := Wi^,k\ (2.32) 

and the abbreviations a U r and r \ a for elements of Q, defined by 

(crUr)^,fe := a^_fcUr^,fc, a = ^\T^ aUr = ^ (2.33) 

If (7 G r^s and r G r2u, we also need ar G f^s+u, defined by 

i(yr)^^k ■= (yuj,k U (s + r^,fe) where g G s + r^,fc ^ g - s G r^,fc (2.34) 

Remark that a function g on fi^ is naturally made into a function on fig+u 
by, using the notation ()2.34|) . 

g{(7T) := g{a), cr G fi^, r G fi„ (2.35) 

13 



2.2.2 Constructing a pathspace measure 

Write the weak coupling generator ()2.7|1 as 

£ = /:o + $^X,fe (2.36) 

LU,k 

with 

J.A-)--=Vkiuj)V:^k-V^,k (2.37) 

and 

LUyk 

Consider Wticr) : B{H) -^ BiTi) as the completely positive map depending 
on cr G fi 

Wt{a) := /n,(a)e*^^°Xi,fc,e(*^~*^)^« . . . e(*i^r*i^i-i)^«X|,|,fc|.,e(*-*i'^i)^» (2.38) 

with Jqj the indicator function of fij C fi and with the indices {u!i,ki),i = 
1, . . . , |cr| referring to the representation ()2.28|) of a. To verify the complete 
positivity of ()2.38|) . rewrite Cq as 

£o(-) = S- + -S\ S = iHj~ lj2'^kiuj)V:^kV.,k (2.39) 

u>,k 

which yields, 

e*^«(.) = e*^-(e*^)* (2.40) 

Complete positivity of ()2.37|) is obvious from its definition. The Dyson ex- 
pansion of e*^, corresponding to the splitting ()2.36p . reads 

e*^= f daWtia) (2.41) 

Jut 

That expression induces a 'path space measure', or a notion of 'quantum 
trajectories' on fl. 

Lemma 2.7. Choose fi G TiTi). Let E G Qt be measurable and define 

P^,,(E):= [ dcrTT[fiWt{a)l] (2.42) 

Je 

Then (P^,t)igK+ o,re a consistent family of probability measures on {Qt)teR+) 
i.e. for a measurable function g on Qt, 

I d^^^t{<y)g{<y) = f dF^Aa)g{a), s > t (2.43) 

where g is extended to Vts as in ()2.35j) . 
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Thus we obtain a new probability measure P^ on Q by the Kolmogorov 
extension theorem, for t > and a function g on fit, 

[ dF^{a)gia) = [ dF^,t{^)g{a) (2.44) 

Jn Jut 

where we used again the extension as in ()2.35|) . The expectation with respect 
to these measures is denoted 

E^,i on fii, E^ on ft (2.45) 

These probabihty measures are often called 'quantum counting processes', 
see |BGMn 4, BMK03]. 

2.3 Results 

We define the integrated entropy current w^ up to time t as a function on Q: 



w" 



j.k 



with n^ f, as in ()2.H2j) . In what follows, we denote by p the stationary state 
for e^^, which is unique by Assumption Al. For k G C, we write 

e{n) := lim - logEje"'^"'*] (2.47) 

it+oo t 

if it exists. Then, e{K) of course depends on all model parameters, i.e. on 
Hs,Vk,r]k. We introduce e^(K) which is derived from the model with new 
parameters Hg, Vi^,r]l = rjk, see ()2.12j) . 

Now we can already formulate the main result of the paper: 

Proposition 2.8. Fluctuation Theorem 

Assume Al. Let w* he defined by \2.4(^ - There is an open set U E C 
containing the real line, M C W, such that for all k G W, the limit 

ein) := lim ^ logEp[e-"'"'] (2.48) 

tt+oo t 

exists and the function k -^ e(/t) is analytic on lA. Moreover, 

e{K) = e\l - k) (2.49) 

// also A3 holds, then e{K,) = e^{K) and 

e{K) = e(l - k) (2.50) 
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We list some consequences of the fluctuation relations ()2.49|) and ()2.50|) . 

Proposition 2.9. Strict Positivity of the Entropy Production 

Assume Al, then 

A2 holds ^ lim -Ep[u;*] > (2.51) 

For the next proposition we introduce energy functions n\ on fi; 

n\:=-Y,^nlk (2.52) 

uieF 

Proposition 2.10. Green-Kubo Relations 

Assume Al and fix some /3 > 0. Let for k, k' G K: 

Lk,k'{P) := iJ- lim ^£^[4] (2.53) 



d(3k' *T+oo t 



f3i=...=f3^jf^=f3 



and similarly the time-reversed coefficient L^. ^,, obtained by starting with H^ 
and Vj! . Then, 

L,AP) + Ll,,{(3) = /? lim iE,K4'] (2.54) 

tj+oo t 

If also A3 holds, then 

Lk,k' = l^hm lEp[nini,] (2.55) 

with Onsager reciprocity 

Lk,k' = Lk',k (2.56) 

2.4 The quantum model: a dilation of the semigroup 

2.4.1 Heuristics 

In the next section we construct a unitary evolution, which is our basic quan- 
tum model. This type of unitary evolutions is generally known as solutions 
of quantum stochastic differential equations, introduced in |HP84j . 



For the readers who are familiar with stochastic calculus, we briefly state 
how our evolution would look in traditional notation. Recommended refer- 
ences are |Par92[ lAttO.Sj for quantum stochastic calculus and |Dern3j for the 
formalism of second quantization. 
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For all tu G F and keK, let (L2(M+))^ ^ be a copy of L^ 



We con- 



sider the bosonic Fock space (F^ denotes symmetrized second quantization) 



7^ = F., 



Lj,k 



(L^(M+)) 



uj,k 



u),k 



(i=(«")).^. 



(2.57) 



and think of dA* ^ ^ with t G M^ as the creation operator on F^ (L^(R+))^ ^ 

creating the "wavef unction" X\t,t+dt\ (the indicator function of the interva 
[t, t+dt]). We now write a Quantum Stochastic Differential Equation (QSDE) 

on B{n®n): 

dU, = 5^r^f(a.)(K.,fcdA;fc,-K:,fcdA.,M)Ut (2.58) 

uj,k 

- (if/;dt-^J]r7fc(^)K:,fcK.,fcdt)Ut, Uo = l®l 
\ i^fc / 

Of course, the intuitive definitions given here, do not suffice to give mean- 
ing to this expression. We content ourselves with stating that ()2.58|) defines 
a unitary evolution Ut, which we will now rigourously construct by using 
Maassen's approach of integral kernels [Maa84j . 



2.4.2 Construction of the unitary evolution U^ 

Recall the Guichardet spaces Vtt and f2, introduced in section 
for (cr, r) G f2 X fi the ordered sequence of times (ti, . . . , t„) as 



and define 
\t\ (2.59) 



{ti, . . . , t„} = U^,fc(o-^,/fc U r^,fc) and < ti < . . . < t„ n = \a 

We define the integral kernel m^ : i7 x f] ^^ B{7i): 

ut{a, t) = In^.nA^, r)e(*-*")^Z„e(*"-*"-)^Z„_i . . . Z^e^'^-'^^'' Z,e'^'' (2.60) 

with IntxUt the indicator function of Qt x ^t, S G B{T-C) as in p.39|l and for 
3 



l,...,n 



if tj G cr^,fc 
if tj G r^,fc 



(2.61) 



1/2/ \T/ 

Finally, let 

J^:=L'^{n,n,da)^n®L'^in,da) (2.62) 

Remark that in a natural way, we have JF ^ 7i (g) 7?. with 7?. as defined in 
(jTKTfl . Take / G L2(fi,H,da) and define 



(u./)(0 = E/«*(^'^)/((^\^) 



Ur)c?r 



(2.63) 



17 



In j Maa84j . one proves that this \Jt is unitary and that it solves the QSDE 

The unitary family Uj, thus defined, is not a group, but a so-called cocy- 
cle; physically this corresponds to an interaction picture and it can be made 
into a group by multiplying it with a well chosen 'free evolution'. Note that 
by taking each V^ ^ = or V^ = in ()2.2p the subsystem decouples from the 
reservoir and ()2.63|) reduces to 

Ui = 1 ® 1 (2.64) 

This follows since the kernel Ut{cr, r) in ()2.60|) vanishes except for cr = r = 
and S reduces to 0. 

Remark that in ()2.58|1 or ()2.(JHj) . the reservoirs are now not only labeled by 
k E K, as in the original physical picture, but also by u; G -F; each transition 
has its own mathematical reservoir. To formulate our results, we also need to 
specify the state. Define the one-dimensional vacuum projection 1^ G B{J-') 

(■\ f\( \ i f(^) when o- = . . 

^^^^)(^) = \0 whena^0 ^^.GS) 

Our reference state is 

p®l0on7^®J^ (2.66) 

where p is the unique stationary state of e^^, see Remark [2. 11 

Note that the state p ® I0 is not invariant under the dynamics, only its 

restriction to 7i is invariant (see also ()2.74j) ). Hence, technically, it is quite 

different from a NESS as in p.2j) . 

We will abbreviate the Heisenberg dynamics as 

it(G) := U,*GUt, GeBin® T) (2.67) 

with Ut as in (J21S1). Let for each A; G K, t > 0, N*^ G B{T) be the energy 
operators 

(Ni/)(a)=4(a)/(a) (2.68) 

with n\ as defined in ()2.52|) . 

We also define a quantity which we interpret as the total energy of sub- 
system plus reservoirs 

N*:=i7s + Xl^fc' ^^^B{n®T) (2.69) 

This interpretation is backed by Proposition 12.131 
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These 'energies' should be understood as renormahzed quantities, of which 
the (infinite) equihbrium energy of the reservoirs was subtracted. This in- 
terpretation is confirmed by remarking that at time s = 0, these 'energies' 
equal 0: for all continuous functions g, 

Tr [l0(7a=oNl)] = Tr [l^giK)] = ^(0) for all fc G iT, t > (2.70) 

2.4.3 Connection of the QSDE with the counting process 

The connection of the QSDE with the 'quantum trajectories' is provided by 
the following lemma, which we will not prove. It can be found for example in 
|BMK03l IBGM04J and it is easy to derive starting from ()2.63|) and remarking 
that 

>Vt(a)(-)=<(a,0)-Ui(^,0) (2.71) 

Lemma 2.11. Let E G fl be measurable (as for \2.S0\) . Denote by 1e the 
ortogonal projection 

1e : L'^i^) -^ L^{,E) (2.72) 

and recall I0 from fl2.65p . Then, for all A G -B(7i) 



Tr^[l0ji(A®l£;)]= / daWt{a)A 

'E 



(2.73) 



where Trjr denotes the partial trace over T . 

The formula ()2.(i3j) actually defines a dilation of the semigroup e*"^. To 
see this, take E = ^^ then ()2.73|1 reads 



Tr^[l0j,(A®l)] 



d(T>Vi(a)A = e*^v4 



(2.74) 



Another useful consequence of lemma 12.111 is the connection between the 
energy operators in ()2.68p and the functions ()2.52|) . 

Proposition 2.12. Let ki, . . . , ki, ti, . . . ,te and gi, . . . ,g£ be finite (i < 00^ 
sequences of, respectively, elements of K , M"*" and continuous functions, and 
let fi G T{7i), then 



Tr 



{^^®M)X{9i{Ju{^l)) 



i=l 



E,, 



9^(11^^ 



i=l 



(2.75) 



Again, we do not give a complete proof and we refer to |BGM04( IBMK03J . 
Proposition 12.121 follows from lemma 12.111 by using that for all t > s and 
keK 

j,{Nl)=UNl) (2.76) 

and that the family {N^ t > 0, A; G K} is commutative. 
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2.5 Results within the quantum picture 

First, we show that the energy (see 12.691 is conserved. 

Proposition 2.13. Let N* be as in \2.69\) . For all continuous functions g: 
Tr [{p ® l0)^(j,(N*))] = Tr [{p Wg^] = Tr [pgiHs)] (2.77) 

The change of entropy in the environment up to time t is 

W*:=^/3,j,(Nt) (2.78) 

k&K 

and its 'steady state expectation' is the entropy production. Our main result 
is a fluctuation theorem for W*. 

Proposition 2.14. Assume Al. Let W* he defined as in \2.7t^ . There is 
an open setUEC containing the real line, M C U, such that for all k ElA, 
the limit 

eU) := lim - logTr \ip ® l0)e-''^*l (2.79) 

tt+oo t L J 

exists and the function n —>■ e{K,) is analytic on U. Let e{K) by defined as in 
!24^ . Then, 

e{K) = e{K) (2.80) 

on U and thus all statements in Proposition \2. SI carry over to c^k). 

From p ® l0 we deduce probability measures Tj on M. Let A C M be mea- 
surable, then 

Ti(A) = Tr[(p®l0)lA(W*)] (2.81) 

where lyi(W*) is the spectral projection on A associated to W*. Via Legendre- 
transformation ()2.50|) implies 

1 , dTA—a) 
- lim - log *\ / = a (2.82) 

which is (jl.l|) . 

In the same way as in Proposition I2.14[ Propositions ()2.10|1 and ()2.9|1 
carry over the quantum picture; for concreteness we give the analogue of 
Proposition (jTTT)|l . 

Proposition 2.15. Assume Al and fix some /9 > 0. Let for k, k' G K: 

Lk,k'{/3) := ^ lim J Tr [{p ® U)Ni)] (2.83) 
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and similarly the time-reversed coefficient Lf. ^,, obtained by starting with if| 



andV^. 



L,,,i/3) + Ll,m = /5 lim - Tr [{p ® M)KK] (2-84) 

tj+oo Z 

If also A3 holds, then 

h,k' = Ip lim ^ Tr [{p ® l0)NiN* ,] (2.85) 

with Onsager reciprocity 

Lk,k' = Lk'± (2.86) 

Recall Lkk> from Proposition \2.1l\ Then for all k & K, 

Lk,k> = Lk^k' (2.87) 

Remark that Propositions 12. 14l and l2?T5l follow immediately from Propo- 
sitions 12.81 and 12.101 by application of Proposition 12.121 



3 Discussion 

3.1 Entropy production for Markov processes 

It is well known that the weak coupling generator is 'classical' in the sense 
that the commutant algebra Ad '■= {A G B{H) | [^, i?s] = O} is invariant. 
In case the Hamiltonian Hs is non-degenerate and only then, Ad is a com- 
mutative algebra. Then we can construct a Markov process with state space 
A which is the restriction of (the dual of) the semigroup e*^ to Ad — C(A). 
Loosely speaking, let p be the stationary state, f2* := A'^'*' the pathspace up 
to time t, and "P* the pathspace measure (starting from p) of this Markov 
process. The time reversal operation O acts on fi* as {Q^){u) = ^{t — u) for 
^ G fi* and < n < t. 

For such Markov processes describing a nonequilibrium dynamics, we 
dispose of a general strategy for identifying the entropy production. It turns 
out in a lot of interesting cases |Ma.en4l IMNVORI IMae99j that 

log^(0 = ^t(0 + C?(l) (3.1) 

where St{^) is the random variable that one physically identifies as the en- 
tropy production. The second term in the righthand side is non-extensive in 
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time. The algorithm allows to derive (jl.lj) from ()3.1|) . 

Since we also have a Markov generator, we can apply the same scheme to 
our setup. ^ To evaluate the result, we however need a physical notion of en- 
tropy production in our model. As mentioned earlier, such a notion is rather 
unambiguous here, see also |LS78| I.TP02bj : 

(current into /c'th reservoir) = Tr [pCkHs\ (3.2) 

But the mean entropy production based on these currents is not equal to 
the expectation value of ()3.1|) : 



V/3fcTr[pAi7s]^lim^Ept 

■'— ' tToo t 

keK 






dQVj, 



(3.3) 



For example, take two reservoirs {k = L(left), i?(right)) and let Refl : BiH) — > 
BiTi) stand for the involution which models left-right reflection. Assume that 
Hs is non-degenerate and that for all a; G M, 

ReRHs = Hs, ReWL = VR, e^^"/V(a;) = e^«"/V(a;) (3-4) 

Hence all parameters are left-right symmetric, except the inverse tempera- 
tures Pl, Pr- (Actually, these assumptions are inconsistent; if Hs is left-right 
symmetric, then it must be degenerate. However, one can introduce an ar- 
bitrarily small symmetry breaking which will generically lift the degeneracy, 
such that the our reasoning still applies.) One checks that 



1 P rr-\ 1 eoGAoCsplfs , V^ , 



em&AmCspHs 

m 

= 0{l)-{(3L + (3R)J2u;i (3.5) 

j=i 

where the sets Aq, A^ and the sequence u^-^ of energy jumps are derived 
from ^, and moreover | Y^^i^il — ll-^sll- This means that in this particu- 
lar left-right symmetric case, ()3.5p is bounded, independently of t for every 
^, and hence the RHS of ()3.3j) vanishes, which disqualifies it as "entropy 
production" . 

^Very recently, a paper IEMQ6J appeared where exactly this is done: one derives a 
fluctuation theorem for log ^jspr as in (|3.1(l . Since the authors consider mainly examples 
involving one reservoir, they do not run into the difficulty described here. 



22 



This trivial remark shows that it is not enough to look at the semigroup 
e^'-' to identify the entropy production. Instead we use more input; we cer- 
tainly use the fact that C = J2keK ^k where the index k runs over the 
different reservoirs but moreover, with the unraveling of the generator, Sec- 
tion I2.2| comes an intuitive interpretation of the various terms. That can 
be contrasted with results by V. Jaksic and C. A. Fillet, where one actually 
proves that quantities like Tr [pCkH^ , cfr. ()3.2|) are limits of currents in the 
original microscopic Hamiltonian model. Of course, we take care that our 
choices are consistent with that result. However, for the higher-order fluctu- 
ations, we do not know; we just make a choice which looks very natural. At 
present, we do not give arguments that for a class of reasonable functions g 

Ep [g{w')\ (3.6) 

is indeed the limit of some fluctuation of dissipated heat in the microscopic 
model. (Although |DR06j points in that direction, see also point 2 in section 

USD 

Another choice for the higher order fluctuations is discussed in Section 1^121 
It is exactly here that lies the role of the dilation with quantum stochastic 
evolutions. If one takes that quantum model as a starting point, then one can 
derive that ()3.6p is a fluctuation of the dissipated heat. To our knowledge, 
that is the only quantum model in which one can study the fluctuations of 
the dissipated heat. 

On the other hand, one can also make a classical dilation of the semi- 
group and in fact, this is exactly what we do in Section 12.21 Yet, there is 
a technical difference between the cases of degenerate and non-degenerate 
system Hamiltonians iJs- If le(-f^s) is one-dimensional for e G spi^s, and in 
addition, for a nonzero cu G F, e is the unique element of spiles such that 
e — a; G spi^s, then we have the following form of Markovianness: If a cr G fi 
contains a;, i.e. 

a = (ToTCTi, (To G Qto, <^i £ ^ti7 T = (^0, ^5 k) for some k E K, to, ti > 

(3.7) 
then 

dP,(a) = Tr[pWt,+tA^)l] 

= Tr [pWto(aor)l] daodr x Tr [U{Hs)Wt,{ai)l] dai 

= dPp(aor)dPi,(Hs)((Ti) (3.8) 

In words, a one-dimensional spectral subspace erases memory. That does not 
work in the degenerate case. 
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3.2 Integrated currents within the semigroup approach 

Starting from ()3.2|) one could define the integrated currents Jk,t £ BiTi) as 

\t= I dMe"^(A(ifs)) (3.9) 

and study tlieir fluctuations. One can ask whether these fluctuations coincide 
with these in our model? The answer is partially positive because 

Proposition 3.1. Take for all k & K : j3k = {3 for a certain f3 and let pp 
be the stationary state for e*^ as in Remark \2.'A For all k,k' & K and all 
u>0, 

d^^^»^^%=o..=u = -T^ [pisC,iHs)e-'{Ck{Hs))] (3.10) 

which gives a relation between Proposition 12 . 1 01 and the Green-Kubo relation 
in |LS99j . Also the averages coincide, leading to ()3.2|) . However it is not true 
that for a reasonable class of functions g 

E,[(7(nt)] = Tr [p(7(4,t)] (3.11) 

So the mean entropy production and the Green-Kubo-formula can correctly 
be expressed in terms of the operators Jk,t, but higher order fluctuations of 
the dissipated heat cannot. 

3.3 Connection to microscopic dynamics 

We know of three derivations in the literature of the stochastic evolution 
()2.(i3|l or ()2.58j) from a microscopic setup: 

1 Stochastic Limit 

Accardi et al. prove in |AFL90j that the weak coupling limit can be extended 
to the total evolution of subsystem observables. Let U^ be the evolution (in 
the interaction picture) on the total system with A the coupling between 
subsystem and reservoirs. Then, in a certain sense, 

U\/,.iS®in),. -^V;iS®l)Vt (3.12) 

whereas the traditional weak coupling limit only speaks about convergence 
in expectation of the left-hand side. The unitary Uf is the solution of ()2.58j) . 
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2 Stochastic Limit Revisited 

In [DR06J, the approach of |AFL90j (mentioned above) was simplified. By in- 



troducing a unitary map Jx acting on the reservoirs, we get for all continuous 
functions g 

s - lim JlUU.Jx = Ui, s - \imJlg{H^,)Jx = giN,) (3.13) 

where s — lim denotes strong operator convergence and H^. is the generator 
of the dynamics in the uncoupled fc'th reservoir. This suggests that one can 
study the fiuctuations of the reservoir energies by looking at the number op- 
erators Nfc in the model reservoirs, exactly as we do in the present paper. 

3 Repeated Interactions 

In |AP03j . Attal and Pautrat describe a subsystem with Hilbertspace TC 
interacting repeatedly for a time h with a small reservoir with Hilbertspace 
TZ. After each time h, TZ is replaced by an identical copy. This procedure 
ensures that at any time, the subsystem sees a 'fresh' reservoir. In the limit 
h —>■ the dynamics (in the interaction picture) converges in a certain sense 
to the solution of a QSDE. One can choose a particular QSDE by tuning the 
parameters of the interaction. Assume that 

7^ = ®c.,fc7^^,fc (3.14) 

Each TZ^^k is 2-dimensional with basevectors {0,uj). Define a^^k on TZcu,k by 

a^,k{^) = 9, a^^kd = (3.15) 

Choose the dynamics on 7i ® 7^ as e~^^^^^^ ioi < t < h with 

H{h) = Hj + ^Yl i^<^.'^<k + V:^,a.,k) (3.16) 

ui,k 

Then, through the limiting procedure of |AP03j . equation ()2.58|) obtains. 



4 Proofs 

4.1 Proof of Lemma 12.71 

From 

d(TWi(a)l = e*^l = 1, W e n : Wt{a)l > (4.1) 
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for all t > 0, it follows that (P^,t)teR+ is indeed a family of probability 
measures for all /i G T(7i). Further, for s,u>0, we have 

Wt{c7)Wt{r) = Wtiar), a e fi„ r G (^, (4.2) 

Together with ()4.1|) . this yields consistency of the family (P^,t)tGM+- 

4.2 Proof of Proposition 12.81 

Define for k G CI^' and t > 0, 

Mi,^ -M^C, Mt;^{a) = J2 ^Aniia) (4.3) 

Our results rely on the following lemma 

Lemma 4.1. Assume Al and let ii G T{Ti,). There is an open setlA C C''^', 
with RI-^I C U such that 

e{-K) := hm -logE„ [e*^*'"] (4.4) 

it+oo t 

is an analytic function on U which does not depend on fi. Moreover, for any 
sequence ki, . . . ,k£ E K, 

hm A . . . ^1 log E, [e^'-j = A . . . ^e(^) (4.5) 

uniformly on compacts. 

Proof. We apply the generalized Perron- Frobenius Theorem 14.31 of the Ap- 
pendix with 

A= /"d(TWt,(a)e^*-^(") (4.6) 

Jn 

for well chosen t^ and r G M)^^. Since for f G M'^', Mt^^f is a real function, the 

map A is completely positive as a linear combination of completely positive 

maps with positive coefficients. 

Below we choose t^ so as to satisfy the non- degeneracy requirement IA-21 of 

the Appendix. 

By faithfulness of the stationary state p, 

e := inf Tr [pP] > (4.7) 

0<PGB(W),P2=P 
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Since the semigroup is ergodic, it follows that there is t^ such that for all 
t>t„ 



sup \\p - e*^>|| < 



(4i 



Since ||£||, ||£o|| < +00, with || • || being the operator norm in B{B(Ti.)), the 
Dyson expansion ()2.41|) is absolutely convergent. Hence, we can find ra G N 
such that 



daW. fa)-e*=^|| < 



a\<n 



3dim?i 



Let m := inf|o-|<n^te,r(cr). For each r G M)^^, decompose 



Mf. 



daWt^Me'""'^ 



e- / daWt^a) 

\a\<n 



(4.9) 



(4.10) 



+ /" da Wt, (a) (e*^'-^(") - e™) + f da W*, (a)e*^*-^(") 

^|o-|<?i J\(T\>n 

and for each pair of non-zero projections P 7^ 0, P' 7^ G B{7i), we have 



Tr 



P / W, (a)e^^*-^P' 



> e" / d(TTr[PWt,(a)P'] 

(T|<n 



>e'"(Tr[Pe*^^P'] - -) 
>e-(Tr[P'p]-^-|)>e-^ 



(4.11) 



This shows that one can apply Theorem 14.31 with A as in ()4.6|1 . Call the 
dominant eigenvalue of A, A(r, t^) and the corresponding strictly positive 
eigenvector v{r). Remark that for each k G C'^' and t G M"*", 



/ 

Jn 



d(T>V,fa)e^*-(") =e*^^ 



where 



u),k 



C^{-) = Coi-) + > :vk{io)e-^'^'-v:, ■ V^,, 



(4.12) 



(4.13) 



This follows by comparing the Dyson expansions (in the same sense as for 
p. 4111 ) corresponding to the left-hand and the right-hand side of ()4.12j) . As 
a consequence, for all r G M'^', Cr has a non-degenerate maximal eigen- 
value A(f) = j-\iaX(f,te) corresponding to the eigenvector v{f). Since v(f) is 
strictly positive, we have Tr [v(f)] > 0, and, for any /i G T(7i), Tr [f (r)/i] > 0. 
This implies 

lim -logTr [/ie^^^ll = A(r) (4.14) 

it+00 t 
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and hence, again by ()4.12|) 

e(r) = A(r) (4.15) 

Since for all Tt G C''^', C^ depends analytically on Tt, perturbation theory 
for isolated eigenvalues gives us for all f G M''^' an open set Ur 3r such that 
for all 'K ^ Ur '■ 

1. There is a unique X(li) G spCf^ such that 

inf{3fJA(7€) - bl |p G sp/:^\ A(k)} > (4.16) 

2. The eigenvector !>(«;), corresponding to A(fi;) satisfies 

inf (3f?Tr[/if(7t)]) >0 (4.17) 



It follows that (pmjl holds for all k G U U-, 



rei 



K\ 



ri 



e{K) = lim - logTr Ifie^^n] = A(7t) (4.18) 

Summarizing, we have for all r G M'^' and /i G TCH) a family of analytic 
functions 

F{t,K) := -logTr [/ie*^^l] (4.19) 

converging pointwise in Ur to the function e(fi;) as t t +oo. 
We recall Montel's Theorem, see e.g. p. 153 of jCon78j : 

Theorem 4.2. Let G G C be open and let (/„)„eN be (i sequence of analytic 
functions G \^ C, then (/„)„eN contains a uniformly convergent on compacts 
subsequence iff the set (/n)neN is locally bounded, i.e., that for each z E G 
there is a r > and M > 0, such that 

\z' -z\<r^yneN: \fn{z')\ < M (4.20) 

For all f G M'^', the family {F{t, 'K))t>to is locally bounded on Ur for large 
enough t^ > 0. This follows from analyticity of C-^ and from the condition 
()4.17|) . Consequently, one can apply Theorem 14.21 for each component of k 
separately. A standard result, e.g. Theorem 2.1, p. 151 in |Con78j states that 
the uniform limit of a sequence of analytic functions is analytic and that all 
derivatives converge. Since this generalizes to the multidimensional variable 
K, e.g. by Hartog's theorem. Lemma f4. II is proven. D 
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Referring again to the representation ()2.28|) . we introduce for a G f2 the 
factor 



ri{a) :-- 



VkA^i) 



(4.21) 



i=l 



Recall the definition of S in ()2.40j) . introduce the time- reversed maps Cq and, 
for t > 0, yyf (i.e., these maps are derived from if| and Vj^) and remark, 
(see also EiOI)) 



J^Co I 



e'' ■ e*^* 



e*^o( 



tTS*T _ tTST 



(4.22) 



Tr 
Tr 
Tr 



Define the operation 6t on Qt as 

6t{u>i, ki,ti] . . .] Un, kn, tn) '■= (— W„, fc„, t — tn] ■ ■ ■ ] —UJi, ki, t — ti) (4.23) 

Calculate 

r\a)Ti[Wt{a)l] (4.24) 

T/* p{'ti+i-U)SY* T/ , p(U+i-U)S*Y , 

■ . . TK;,^„,,^,TTe(*'+-*')^*TTK;,,,T . . . TV^,* ,,TTe(*»+-*»)^TTV:,*^^,,^^^T . . . 

• • • TK!,^^^,,^^^Te(*'+-*')^^*^TK!,^,,T . . . TV^_.,,,Te(*'+-*')^^^Ty_,^,,,^,T . . . 
= v'\9t(T) Tr [WfiOtcr)!] = r/-i(a)e""(") Tr [W^{9tcr)l] 

In the last equahty the KMS-condition ()2.3p was used. The previous equali- 
ties follow from cyclicity of the trace, TT = 1, V^^ = V^-w,fc and ()4.22|1 . Using 
()4.24|) . we calculate by change of integration variables (putting I :- 

Tin), 



1 



diniH 



Ej [e^*--] = /" da Tr [IWt{a)l] e*'*--^'^) = / da Tr [/Wf (a)l] e-^'^'^e-^*--^") 
Jn Jn 

(4.25) 

Since in the limit t | oo, one can replace the initial state / by p, as in ()4.18|) . 

the formula ()4.25j) yields for all k G W as in Lemma 14.11 

e{K) = e\l - -R) with 1 - k := (1 - Ki, . . . , 1 - k\k\) (4.26) 

Finally, Proposition 12.81 follows from ()4.2(jp by putting for some k G C, 

Ki:=K i = l,...,\K\ (4.27) 

thus obtaining Mt^n = f^w^ 
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4.3 Proof of Proposition 12.91 

The nonnegativity of the entropy production follows from Proposition 12 .81 bv 
Jensen's inequality. To get the strict positivity from A2, we first need to 
introduce more notation. 



Let p be the unique stationary state of e*^. We decompose the states p 
and TpT in one-dimensional unnormalized states as 



P 



Up' 



p-p- 



kMWpl P^>0, 1,3 eD 



(4.28) 



ieD 



where p* can stand for p or TpT and D := {!,..., dimTi}. The decompo- 
sition ()4.28|) differs from the spectral decomposition when p is degenerate. 
Remark that there is an arbitrariness in labeling the unnormalized states, 
as well as a possible arbitrariness stemming from degeneracies in p*. We 
partially fix this arbitrariness by asking that 



T{Tp*T)jT = p* 



(4.29) 



This is always possible because the set T{Tp*T)jT,j G D satisfies all the 
requirements of ()4.28|) as a decomposition of p*. Let Cl^ = Qt >^ D x D ioi a 
t > and define the measure P^ by (letting ^f be a measurable function) : 



/ dPi(a)(7(a) = V / daTr 
Jnt ij Jnt 



P^m{(r) 



(TpT), 

II (TpT), 



(4.30) 



where it is understood that cr G fit and i,j G D. In the rest of this section we 
will use this notation without further comments. Positivity of Pt is obvious 
and normalization follows by 






do"Tr 



P^^^t{(T) 



(TpT), 



E 



II (TpT), II J 
daTT[p,Wt{a)l]= /"dPp,t(a) 



(4.31) 



We call Pf the measure, constructed as above, but with Wf replacing 
Wt. Remark that this is not the measure one would obtain by starting from 
Hg, V^ instead of i^s, V^, because then one would also replace p in ()4.30|) by 
p^, the stationary state of C^. 
Define again the operation 6t on Qt as 



e 



t{(r,t,j, 



ncr,J,t) 



(4.32) 
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where the action of dt on Vtt was defined in fl4.23p . 
Consider the function 



5*: a 



s 



H~a) 



log- 



d¥^{9d) 



a 



We upgrade the function w* on Vtt to a function on Vtt as 



w\a)=w\a), a={a,i,j] 



Our strategy will be to prove f Section 14.3.11) that for some u > 



dP„(a)5"(a) > 



i lu 



and then (Section l4. 3. 2|1 that for allt > 



fit 



I (5* 



^t[a)[S\a)-w\a)) <o 



which will lead to the conclusion that for a certain u G 



n Jn 



\(a)t/7"(a) > 



(4.33) 



(4.34) 



(4.35) 



(4.36) 



(4.37) 



where the first equality is checked by arguing as in ()4.3H) . The converse 
statement is proven in Section 14.3.31 

4.3.1 Positivity of S* 

Looking back at the calculation ()4.24|) . one immediately checks that for t > 
and a E Q, 



Tr 


[\p.r'^''^\{TpT),\\ 


= e"'*^'') Tr 


.\pj\ \KTpT)i\ 


and hei 


ice 







(4.38) 



S\a)=w\a)-\og{\\p,\\) + \og{\\p, 



(y = [(y,hj} 



(4.39) 



Note, using ()4.33p and S^{a) = —S^{9tcr), that S** satisfies an exact fiuctua- 
tion symmetry, for t > and «; G C: 



dPi(a)e-"^'(^) = 
fit J fit 



^^(a)e-(^-'^)^*('^) 



(4.40) 
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Remark that f : M. —>■ W : x —* e ^-|-x — lis positive for all x, increasing for 
X > and decreasing for x < 0. A Chebyshev inequality with 6 > yields 



t/':i\ _ / ^m> /:^\/„-S* , nt 



d¥tia)S\a) = / dPt(a)(e-*^ + S'- l)(a) > ie-' + 5 - im\S'\ > 5) 
'nt Jnt 

(4.41) 

Rephrasing ()2.15j) - ()2.1fi|l . there is for m > 0, a i? C f2„, and one-dimensional 

projection P E BiTi) such that 

f daTr [PW„(a)P] > 0, w"(E) =: w ^ (4.42) 

Je 

For any fc G N, we construct 

Qku D E^ := {ai(r2 ...ot, \ai, . . . ,Ok e E} (4.43) 

where the notation 0"icr2, and consequently also o"iO"2---crfc was defined in 
(Hini. We have 

w\E^) = kw, [ daTr [P>V„((7)P] > (4.44) 

Since p is faithful, there are i,j G {1, . . . , dimTi} such that 

/ daTT[piWku{cr)pj]>0 (4.45) 

Je'' 

Since the function S*" — w^ is bounded uniformly in n G M^ (this follows 
e.g. from ()4.39j) ). one can choose k eN and z, j G {1, . . . , dimTi} such that 

Tr [piWkt{(r)pj] = kw + log \\pi\\ - log ||p,|| > 0, a e E'' (4.46) 

This proves that the last expression in 1)4.411] is not zero (after replacing t by 
ku). Hence, ()4.35p is proven. 

4.3.2 Difference between 5* and w^ 

Calculate for t > 



/ dP,(a)log||p,||= Y, 

•^"* i,j&D 

Tr [pe*^ log TpT] = Tr [p log TpT] 



[ daTi 


V^'^'<:^\ 


loglPil 


iTpT] 




(4.47) 
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and 



[ d¥{a)log\\pi\\ =J2 / daTr 



P^Wt{(T) 



log 1 1 Pi 



IIp.II 

= ^Tr[pie*^l]log||p,|| = Tr[plogp] (4.48) 

ieD 

where we used pe*"^ = p and e*^l = 1. 
Hence one gets 

/ dPt(a) ( log((TpT),) - log(p,)) = Tr [p(logp - log TpT)] < (4.49) 

where the last inequahty follows from the nonnegativity of the relative en- 
tropy. 

4.3.3 Strict positivity implies Assumption A2 

We prove that A2 is a necessary condition for a nonzero entropy production. 
First remark that 

/ dFp{a)w\(j) (4.50) 

is extensive in t > 0. This follows from translation invariance (in t) of w^ 
and stationarity of Pp. Hence we can fix t > such that 

drp{a)w\(x)\ > 2dim7^max|/3fccj| (4.51) 

Take a & Vtt satisfying Wt{cr) ^ 0. It follows that one can split t = X]i=i ^j 
and 

(T = rarsTi, n E ^t^, i = 1, 2, 3 (4.52) 

(again the notation ()2.34|) was used) such that 

1. There is a one-dimensional projection P such that 

Tr[PWi,(r2)P] >0 (4.53) 

2. 

|ri| < dim?^, Irg] < dmxiH (4.54) 

Assume that A2 does not hold. It follows that w^'^{t2) = 0. 
Hence, by (ji3i|l . 

\w\a)\ = \w^'{ti)+w^'{t3)\ < 2dimnmax\pk^\ (4.55) 

which is in obvious contradiction with 1)4.511) . 
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4.4 Proof of Proposition 12.101 

This proof is by now quite standard, it can be found e.g. in |LS78j . We 
recall from ()4.5j) in Lemma f4. II that we can interchange the limit t t oo and 
differentiation of k h-i> e(7t). By differentiating relation ()4.2fj|l with respect 
to Kk and to (3k' in k = and f3k^K = P, and interchanging limits and 
derivatives, we arrive at the modified Green-Kubo relation: 

Lk,i + Lli = dhm ^Ep[ninl,] (4.56) 

from which the other statements in Proposition I2.1UI easily follow. 

4.5 Proof of proposition 12.131 

Choose V ETi and G L'^iVtt) such that 



H^v = msv \/k e K : N^0 = m^^ ms, nikaK G I^ (4-57) 

By the definition of Ut, 

\Jt{v(^<Pm = Y. [ dTUt{cx,T)v<l){{^\a)UT), eea (4.58) 

Using [S, Hs] = 0, one checks that Uti^cr, t)v either vanishes or 

Hs{ut{(T,T)v -v) = ^u;(|(7a;,fc| - |r^,fc|)(Mt(a, r)t; -t;) (4.59) 

u),k 

By fg37ll . it follows that 0((^ \ a) U r) = in (H3H1) . unless for dX\ k e K 

^^{\ii^A - |o-c^,/c| + ka;,fc|) = mk (4.60) 

Together with ()4.59j) . this implies 

N*Ut(t; (g) 0) = N*(w (g) 0) = [ms + J^ mk){v ® 0) (4.61) 

k&K 

Since the operators i/s, ^\^k mutually commute, vectors like t> as in 
()4.57|1 furnish a complete set of eigenvectors. This proves the proposition. 
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4.6 Proof of Proposition 13.11 

By expanding the left-hand side of ()3.10|1 in a Dyson expansion, as in ()2.41|) . 
one can evaluate the derivatives, leading to 

■d^^P^«%.=.,..=u = E ^1^2 Tr [pX„.,e"^X„..(l)] (4.62) 

Putting p = pf^, yields ppV^^^k = V^^kPp^'^^"^ ■ Now ()3.10|) follows after some 
reshuffling, using 

K:,fc = V.^,k, r/fc(^) = e-^'=-r/fc(-^), ^a;r7,(^)K:,fcK.,fc = A(/^s) 

(4.63) 

APPENDIX 

Let A be the matrix algebra M„(C) for some n G N, and denote by A^ its 
positive cone, i.e, 

A^ = {x*x \x e A} (A-l) 

An element x G A'^ is called strictly positive (notation: a; > 0) if it is 
invertible. 

Theorem 4.3. Let A : A ^ A be a completely positive linear map, satisfying 

Tr[xAy]>0, x,y e A'^,x y^ 0,y y^ (A-2) 

Then, A has a positive eigenvalue X, such that if p is another eigenvalue, then 
\p\ < A. The eigenvector v ^ A corresponding to A can be chosen strictly 
positive. The eigenvalue A is simple, i.e., as a root of the characteristic 
equation of A it has multiplicity 1. 

The theorem was proven almost in the above form in |EHK78j , (see The- 
orem 4.2 therein). We state (a simplified version of) that theorem and we 
show that the above statement follows from it. We call a positive map on 
A irreducible if 

Vx ^ 0, y ^ G ^+, 3A; G N : Tr [xcp^] > (A-3) 

Theorem 4.4. Let cj) be a positive map such that 

1. (j) preserves the unit 1 G A: 0(1) = 1 
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2. (p satisfies the two-positivity inequality: 

(f){x*x) > (f){x)*(f){x) for all X e A (A-4) 

3. For all k = 1,2,..., 0'^ is irreducible 

Then, has a positive, simple eigenvalue X, such that if fi is another eigen- 
value, then \n\ < A. The eigenvector v G A corresponding to A can be chosen 
strictly positive. 

Another theorem in |EHK78j is (Theorem 2.4, combined with the sen- 
tences following it): 

Theorem 4.5. Let be an irreducible positive linear map on A and let r be 
the spectral radius 

r := sup{|c| I c G sp0} (A-5) 

then there is a unique eigenvector v G A'^ with eigenvalue r. 

To prove Theorem 14. H| we remark that A has the same spectral properties 
as a well-chosen map that satisfies the conditions of Theorem 14.41 
Since A is irreducible, one can apply Theorem 14.51 to find an eigenvector v. 
Because of IA-21 we conclude that f > 0. Let now the map be defined as 

0(x) = -v-'^^A{v^/^xv^/^)v-^/^, xeA (A-6) 

It is clear that 

1 . is completely positive and still satisfies IA-21 

2. 0(1) = 1 

3. sp0 = -spA and also the multiplicities of the eigenvalues are equal 

Hence satisfies the conditions of Theorem 14. 4t since unity-preserving 
completely positive maps satisfy the two-positivity inequality ()A-4jl . Theo- 
rem 14.31 follows. 
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